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INTRODUCTION: MOTIVATION 

In the framework of 'instantaneous approximations' to the Bethe-Salpeter formalism for 
the description of bound states within quantum field theories, depending on the Lorentz 
structure of the Bethe-Salpeter interaction kernel solutions of the (full) Salpeter equation 
with some confining interactions may exhibit certain instabilities [1], possibly related to 
the Klein paradox, signalling the decay of states assumed to be bound by these confining 
interactions, and observed in numerical (variational) studies [1] of the Salpeter equation. 

The (presumably) simplest scenario allowing for the fully analytic investigation of this 
problem is set by the reduced Salpeter equation [2] with harmonic-oscillator interaction. 
In this case, the integral equation of Salpeter simplifies to either an algebraic relation or a 
second-order homogeneous linear ordinary differential equation, immediately accessible 
to standard techniques. There one can hope to be able to decide unambiguously whether 
this setting poses a well-defined (eigenvalue) problem the solutions of which correspond 
to stable bound states associated to real energy eigenvalues that are bounded from below. 

REDUCED SALPETER EQUATION FOR CONFINING 
INTERACTION KERNELS OF HARMONIC-OSCILLATOR TYPE 



Assuming, as usual, the Lorentz structures of the effective couplings of both fermion and 
antifermion to be represented by identical matrices in Dirac space, generically labeled F, 
and denoting the associated Lorentz-scalar interaction function by Vrlp, q) , the so-called 
reduced Salpeter equation [2] describing bound states composed of some fermion and its 
corresponding antifermion (with common mass m and relative momentum p) reads for a 
bound state with mass eigenvalue M in the center-of-momentum frame of the bound state 

(M-2£)4>(p)=A+(p)7o J ^£yr(p,q)r4>(q)rA-(p)ro, (1) 



with the one-particle kinetic energy E and energy projectors A='= (p) defined according to 
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Any solution to this eigenvalue problem is the Salpeter amplitude ^(p) of a bound state. 

Let the Bethe-Salpeter kernel be of convolution type, Vr(p, q) = H^(p — q) ? arising from 
a central potential y(r) , r = |x| , in configuration space. Then, for the harmonic-oscillator 
potential V {r) — ar^ ^ a — a* ^ the reduced Salpeter equation becomes a second-order 
differential equation utilizing the Laplacian acting on states of angular momentum £ = 0, 

d2 2d 
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In order to make contact with related previous analyses [3-6], let us present our line of 
argument for fermion-antifermion bound states of total spin /, parity P = {—iy~^^, and 
charge-conjugation quantum number C = ( — 1)'^, called ^Jj spectroscopically. Due to the 
projectors A"^ (p) on the right-hand side of the reduced Salpeter equation (1), the Salpeter 
amplitudes <I>(p) describing these states contain only one independent component 0(p): 

^>(p) = 20(p)A+(p)r5. 

Somewhat more specifically, we consider pseudoscalar (^Sq) bound states, characterized 
by the spin-parity-charge conjugation assignment = ^. A brief moment of thought 
reveals that the cumbersome instabilities should appear first in the pseudoscalar channel. 
Stripping off all dependence on angular variables [7] converts this "harmonic-oscillator 
reduced Salpeter equation " into the eigenvalue equation of a (Hamiltonian) operator : 

jr0(p)=M0(p). (2) 

It is a straightforward task to work out the explicit form of all the Hamiltonians M' for 
the most popular choices of the Lorentz structure of Bethe-Salpeter kernels (cf. Table 1). 

TABLE 1. Hamiltonian (differential or mere multiplication) operators entering into the eigenvalue 
equation (2) equivalent to the reduced Salpeter equation ( 1 ) with harmonic-oscillator interaction potential 
V{r) =ar^, for several frequently considered Lorentz structures of the Bethe-Salpeter interaction kernel. 

Lorentz structure F^F Jf 

Lorentz scalar 101 \~2E^ E 

n n ^ f2p^ + 3m^ \ 
time-component Lorentz vector y^igjy" 2E+a\ ^^-^ D j 

( n? 1 

Lorentz vector Jix.®'^ 2E + a[—D—-2D 

Lorentz pseudoscalar 75 (g) 75 2E + a ^-^^^ 

Bohm-Joos-Krammer (BJK) [8] 5 (y^t ® 7^ + 75 ® 75 - 1 ® l) 2E-aD 



SPECTRAL PROPERTIES OF THE REDUCED-SALPETER 
HAMILTONIAN OPERATORS 



The spectra of the mukiplication or differential operators ^ for various Dirac structures 
r ® r fixed by the Bethe-S alpeter kernel exhibit the following stability -relevant features : 

• All of our Hamiltonian operators M' are self-adjoint since the differential operators 
D and DE^ as well as the multiplication by any real- valued function define 

self- adjoint operators. Consequently, the entire spectrum of any operator is real. 
For "reasonable" interaction kernels, that is, for kernels only constructed in terms of 
Dirac matrices T subject to Jq Jq = ±r and potential functions VrlP? Q) subject to 
P) = ^r(p, q) , the reality of all eigenvalues M is guaranteed by the relation [7] 
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satisfied by any Salpeter amplitude ^(p) solving the reduced Salpeter equation (1). 
For the Lorentz pseudoscalar, F ® F = 75 C?) 75 , and, if m = 0, for the Lorentz scalar, 
F(8)F = 101, our Hamiltonians form pure multiplication operators, with purely 
continuous spectrum. Bound states do not exist, stability questions thus do not arise. 
For the time-component Lorentz vector, F F = 7^ ® 7^, for the Lorentz structure 
introduced in [8], F(g)F= j (7^ (g) 7^ -1-75(8)75 — 1 1), and, if m = 0, for the Lorentz 
vector, F (g) F = 7^ (g) 7^ , our Hamiltonians form (i = 0) Schrodinger operators 
with a positive, infinitely rising potential y(p) 0° for p ^0°, provided, of course, 
the sign of the coupling a is chosen appropriately. The latter operators have entirely 
discrete spectra bounded from below; all bound states may be expected to be stable. 
Figure 1 illustrates the typical qualitative behaviour of the effective potentials V{p). 




FIGURE 1. Effective potential y (79) in the Hamiltonian for the harmonic-oscillator reduced Salpeter 
equation of time-component Lorentz-vector kernel F ® F = 7" (g) y" , for m — 1 and a = 1 [arbitrary units] . 
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FIGURE 2. Lowest-lying eigenvalues e(M) of the auxiliary Hamiltonians ^ = —D + f/(/7;M), derived 
from the reduced Salpeter equation (1) describing pseudoscalar (^So) fermion-antifermion bound states of 
constituents with mass m=l, which experience a harmonic-oscillator interaction [V{r)^ar^,a — a* ^ 0] 
(a) of the Lorentz-scalar structure FigjF = 1 1 (with the "binding" coupling constant a = — 10) and (b) of 
the Lorentz-vector structure F ® F = y^j (g) 7^ (with a "binding" couphng constant a =10) [arbitrary units] . 



• For m 7^ 0, because of the presence of the nasty differential operators m^E ^ DE^ ^ 
the Hamiltonians in (2) corresponding to both Lorentz scalar, F F = 1 ® 1 , and 
Lorentz vector, Y®Y= y^i^J^ , do not constitute (standard) Schrodinger operators. 
In these cases, however, by a suitable redefinition of the radial amplitudes (p) , the 
(radial) differential equations may be transformed to eigenvalue equations of (i = 0) 
Schrodinger operators ^ of the form J^^ = —D + U{p; M) , which involve effective 
potentials U{p;M) depending on the bound-state mass M as parameter. As might be 
guessed from the form of the corresponding Hamiltonian for the Lorentz scalar 
the transformation required here simply reads ^{p) ^ E ^ (p) . For given M, and the 
appropriate sign of a, our effective potentials U{p;M) are bounded from below and 
behave like U {p;M) ^ oo for p ^ oo. Thus, the spectra of both "auxiliary" operators 
must consist entirely of discrete M-dependent eigenvalues e{M) (cf. Fig. 2). The 
derivatives of the latter eigenvalues with respect to M are, for all M, strictly definite. 
The bound-state masses M, defined by the zeros of the eigenvalues e (M) , must then 
be discrete too. Since these eigenvalues e(M) are strictly decreasing functions of M 
a closer inspection establishes the bound-state masses M to be bounded from below. 

In summary, given the semiboundedness of all the Hamiltonians J^f entering in the radial 
equations and having established the discreteness of their spectra for harmonic-oscillator 

couplings of appropriate sign, our harmonic-oscillator reduced Salpeter equation poses, 
at least for a very wide class of Lorentz structures, a well-defined problem, with solutions 
describing stable bound states related to some real discrete spectra bounded from below. 



GENERALIZATION TO THE (FULL) SALPETER EQUATION 



It goes without saying that a similar study may be envisaged for the full Salpeter equation 
(as before, for bound-state constituents of equal mass in the rest frame of the bound state) 



There, however, the spectral analysis will be more complicated for the following reasons: 

• Although the squares of all the mass eigenvalues M are guaranteed to be real [9] 
the spectrum is, in general, not necessarily real and, even in those cases where it can 
be shown to be real, it is certainly not bounded from below [9]. In particular, for the 
perhaps most important example of Bethe-Salpeter kernels involving only potential 
functions Vr(p,q) satisfying Vp(p,q) = Vr(p,q) = Vr(q,p) and coupling matrices 
F satisfying Jq F^ Jq = ±F the spectrum of mass eigenvalues M in the complex-M 
plane consists of real opposite-sign pairs (M, —M) and imaginary points M = —M* . 

• "FuU-Salpeter amplitudes" 4>(p) , that is, solutions of the full Salpeter equation (3), 
have more than one independent components. Thus, any full Salpeter equation with 
harmonic-oscillator interaction translates to a system of more than one second-order 
differential equations or, equivalently, a single differential equation of higher order. 





(3) 



The Salpeter amplitude for states, e.g., involves two independent components, 0i , (fe: 



4>(p) 



(p) ^ ^ + 02 (P) 
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In particular, for interaction kernels of Lorentz- scalar (r(g)r = \®\) or time-component 
Lorentz-vector (F^F = y*' ® structure, the full Salpeter equation (3), after getting rid 
of all angular variables, becomes equivalent to the fourth-order differential equation [10] 

4E^-^[G{p^ + 2m^)D + pDp-2\+^{m^D + opDp-2o) -d\^2{p) 

with a sign factor a (rendering possible a simultaneous study of both structures) given by 

forY®T="P®"P (time-component Lorentz-vector interactions) , 
— 1 forr(8)r=l(8)l (Lorentz-scalar interactions) . 



A notable exception is the Salpeter equation (3) with a harmonic-oscillator interaction of 
the BJK [8] Lorentz structure T®T =^ {y^ ®7'^-|-75(E)75 — l<E)l): here one still arrives 
at a second-order differential equation that can be expressed in two equivalent ways [10], 

A{E^-aDE)(^x{p)=M^(^i{p) or A {E^ - aE D) (h{p) = (h{p) , 

since 2 £ 0i (p) = M 02 (p) . Again, the spectrum of bound-state masses M is discrete [10]. 
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